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Abstract 

The continuum 3d SU(2)xU(l)+Higgs theory is an effective theory for a large class of 
4d high-temperature gauge theories, including the minimal standard model and some 
of its super symmetric extensions. We study the effects of the U(l) subgroup using 
lattice Monte Carlo techniques. When g' 2 /g 2 is increased from the zero corresponding 
to pure SU(2)+Higgs theory, the phase transition gets stronger. However, the increase 
in the strength is close to what is expected perturbatively, and the qualitative features 
of the phase diagram remain the same as for g' 2 = 0. In particular, the first order 
transition still disappears for tuh > ttih,c- We measure the photon mass and mixing 
angle, and find that the mass vanishes in both phases within the statistical errors. 
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1 Introduction 



For reliable computations in the cosmology of the very early universe, it is indispensable 
to know the equation of state of the matter governed by the laws of the one correct - 
so far unknown - physical electroweak theory. For a large class of candidate theories 
an effective theory, the 3d SU(2)xU(l) + fundamental Higgs theory, can be derived by 
analytic perturbative computations |], |2[. This effective theory gives an approximate 
but in many cases accurate description of the electroweak phase transition. The non- 
perturbative effects always encountered in a finite T context are isolated in the effective 
theory, which has to be solved numerically. Up to now, all non-perturbative efforts [3- 
8] were concentrated on an even simpler theory, without the U(l) factor. The 4d 
lattice simulations of the bosonic sector of the standard model || also do not take 
into account the U(l) group. (The present status of lattice Monte Carlo studies of the 
SU(2)+Higgs model both in the finite temperature 4d theory and in the 3d effective 
theory is reviewed in [fL0| .) The purpose of this article is to complete the study by 
including the U(l) subgroup. 

At first sight, the inclusion of the U(l) factor is harmless and may be done pertur- 
batively. However, this need not be so in the vicinity of the phase transition. The 
dimensionless high-temperature expansion parameter associated with the U(l) group 
is ~ g' 2 T/m, where g' is the U(l) gauge coupling and m is a typical mass of a (hy- 
per) charged particle in 3d. This parameter is indeed small deep in the broken phase 
(where m is the W boson mass, m ~ gv) and high in the symmetric phase (where m 
is the mass of a hypercharged scalar excitation, m ~ gT). However, near the critical 
temperature the scalar mass may be small and perturbation theory breaks down. Thus 
a Monte Carlo study is needed. 

The inclusion of the U(l) factor also changes qualitatively the mass spectrum of the 
SU(2)+Higgs theory, as there is now an excitation (the photon) which is perturbatively 
massless in both phases. There are several interesting questions related to this exci- 
tation, for instance whether there might be non-perturbative mass generation for it 
and how the Weinberg mixing angle (defined in a suitable gauge independent manner) 
behaves in the vicinity of the phase transition. 

It is important to keep in mind the separation of the two stages in the effective 
field theory approach. The first stage establishes the values of the parameters of the 
3d effective theory in terms of the physical 4d parameters by means of a perturbative 
computation. This has already been done for the minimal standard model (MSM) in 
|l| and for the MSSM in @j. The second stage is a non-perturbative lattice study of 
the 3d theory as such, performed in the present paper. 

Since the simulations are rather demanding in computer time, we do not chart the 
entire critical surface but perform simulations only at two new points of the parameter 
space. The nonzero value of the U(l) gauge coupling gf = 0.3g 2 comes from the 
ratio of the physical values of m\y, in the minimal standard model. For the scalar 
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self-coupling we choose two values: 

• A3 = 0.06444(^3, which in the notation of || corresponds to m* H = 60 GeV (for 
the actual physical Higgs pole mass ran in different 4d theories, see Sec. [|). This 
value of the "Higgs mass" was the most precisely studied case in the SU(2)+Higgs 
model. The first order nature of the phase transition here is clear. 

• A 3 = 0.62402pg, corresponding to m* H = 180 GeV. In the SU(2)+Higgs model 
first or second order phase transitions are absent then ||. We find that the same 
statement is true when the U(l) interactions are added. 

The paper is organized as follows. Sec. 2 is the formulation of the 3d problem in the 
continuum, Sec. 3 gives examples of 4d— *-3d connections and in Sec. 4 we describe how 
the theory is discretized. Some perturbative results are in Sec. 5. Sec. 6 contains the 
lattice results and a comparison with perturbation theory, and the conclusions are in 
Sec. 7. 



2 Formulation of the problem in 3d 

The Lagrangian of the effective theory (we keep only super-renormalizable interactions; 
the accuracy of this approximation in the MSM is discussed in |iT| ) is: 

L = ±F$F* + Iflyfly + (A0) f D4 + m^V + A 3 (0V) 2 , (2.1) 

where 



Bij = diBj — djBi 

Di = di + ig 3 A, + ig' 3 Bi/2, A { = l -o a A1, 




The SU(2)xU(l) local gauge transformation is 

(j)(x) -> e ia(x) G{x)<f){x). (2.3) 

The four parameters of the theory, g%, g'^, A3 and m\ are definite computable functions 
of the underlying 4d parameters and the temperature. 

From the 3d point of view, the coupling constants g%, gf and A3 are renormalization 
group invariant, whereas the mass parameter 777J gets renormalized at 2-loop level. In 
the MS scheme, the renormalized part of m\ is of the form 

^) = #^log— , (2-4) 
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where 

h m = ~ f<&? ~ + 9A 3 ^ 3 2 + 3A 3 2 3 2 - 12A^, (2.5) 

and A m is a constant specifying the theory. 

All the four parameters of the 3d SU(2)xU(l)+Higgs theory are dimensionful. One 
may measure all the dimensionful observables in terms of one of the coupling constants, 
say g\. The dynamics then depends on the three dimensionless parameters x, y and z, 
defined as 

93 93 93 

These quantities are renormalization group invariant. 

The question we are going to address is: " How does the phase structure and how do 
the gauge invariant operator expectation values depend on z > 0?" When z = the 
U(l) interactions are decoupled and we get an SU(2)+Higgs theory, the phase structure 
of which has been studied in [[| f|. In the opposite limit, z — > oo with x/z fixed, we 
get a 3d U(l)+Higgs model with global SU(2) symmetry. The physical value of z lies 
in between and is related to the Weinberg angle, z ~ m|/m^ -1« 0.3. 

The main characteristic of the 3d SU(2)+Higgs theory is the existence of a critical 
line y = y c {x) for < x < x c ~ 1/8, along which the system has a first order 
transition. The transition gets weaker with increasing x and eventually terminates in a 
second order transition possibly of the Ising universality class at x c . At x > x c there is 
no transition or a transition of higher than second order. If there is no transition, the 
two phases are analytically connected, which is possible since there is no local gauge 
invariant order parameter which would vanish in either of them [12-15]. The signal of 
the first order transition is that along y = y c (x) the system can coexist in two phases 
with different values of various gauge invariant operators like (</>V)- 

Correspondingly, for the SU(2) xU(l)+Higgs theory one expects that there is a 
critical surface y = y c (x,z), which for z = coincides with the previous curve 
y = y c (x). For larger x the surface should terminate in a line of second order transi- 
tions, x = x c (z), x c (0) ~ 1/8 (Fig.l). Indeed, the statement about the non-existence 
of a local gauge invariant order parameter is unchanged. This expectation has to be 
confirmed numerically, and the precise magnitude of the effects of making z > on var- 
ious physical quantities have to be determined. Expressing everything in dimensionless 
form by scaling by powers of gf , the quantities we are interested in are 

• The critical curve y = y c (x, z), 

• The jump A£ 3 of the order parameter like quantity £ 3 = ((p^ <p(g 3 )) / g 3 between the 
broken and symmetric phases at y c . In perturbation theory, A£ 3 ~ \$ ) {y c ) f g\, 
where <pb is the location of the broken minimum in, say, the Landau gauge, 
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Figure 1: The expected schematic phase diagram of the SU(2)xU(l)+Higgs theory in 
the space of the parameters of eq.( |2.6|) . There is a 1st order transition which terminates 
in a line of 2nd order transitions. 

• The interface tension a 3 , defined in perturbation theory by 

^3 = j* b ' 9Z d{<P/g^2V{<P/g,)/gl (2.7) 

A non-perturbative study of the dynamics of the theory is based on a study of the 
expectation values and correlators of low-dimensional gauge invariant operators, so 
for completeness we will review these here. It is convenient to introduce a matrix 
parametrisation of the Higgs field by writing 

Under an SU(2)xU(l) gauge transformation $ transforms according to 

G{x)^e~ ia{ - X > 3 . (2.9) 

The global transformation 

$(ac) -> §{x)C- 1 (2.10) 
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where C is an SU(2) matrix, will be seen to be an isospin custodial rotation of the 
vector states W*. 

In order of increasing dimensionality the list of gauge-invariant operators is: 

• Dim = 1: <fr^<f). Note that (jy<fi = 0; there is only one singlet. 

• Dim = 3/2: By, the U(l) field. 

• Dim = 2: there are two independent gauge invariant operators, namely a vector 
4>^Di(j) — [Di<j))^<j) and a scalar (0V) 2 - F° r the SU(2)+Higgs case (z = 0) there 
are two extra vector operators. In terms of the matrix representation (|2.8|) the 
vector operators invariant under the SU(2) subgroup are 

= Uti [$t - (D^)^]a a . (2.11) 

Under the U(l) gauge transformation ( |2.9|) (which is a rotation around the 3-axis 
in isospace) 

W? ^>iTr&D i $(e- i6 " a (T a e i9 ' T3 ), (2.12) 
so that W}' 2 are not gauge-invariant but Wf is. Note that 

Tr [$ f A$ + (A$) t( £] = 2^(0 f 0) (2.13) 

and does not bring in anything new. We remark also that in the SU(2)xU(l) 
theory the custodial symmetry is explicitly broken. Under the custodial trans- 
formation ( |2.10| ) the operators in (p. 11] ) transform as vectors in the space of 
a = 1, 2, 3, 

W t a -> zTr tfD&C-^aC, (2.14) 
and thus a mixing of gauge-invariant and gauge-noninvariant operators occurs. 

• Dim = 5/2: ^F ij( p, fl^By. 

• Dim = 3: Terms in the action. 

• Dim = 7/2: The parity violating J PC = 0~ + terms 

P 1 = te ijk [^D k( p-(D k( l))U]B tJ , (2.15) 
P* = ie^F^D^-iD^F^}. (2.16) 

The list of parity odd operators can be extended by two non-local operators, 

e i2 ^cs, (2.17) 
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where 

N C S2 = / d 3 xe ijk Tr - -igAiAjA^j (2.18) 

and f 

N C sx = j^J ePx€i jk BiB jk . (2.19) 

In this paper we mainly use the scalar operator the U(l) operator B^ and the 
vector operator Wf. 

3 4d — > 3d relation 

As emphasized earlier, it is important to separate the two stages of solving the finite 
T problem: relating the 4d and 3d theories and discussing the 3d theory as such. The 
mapping here is many — > few; many different 4d theories with different parameter 
values can correspond to a single set of 3d parameters x, z. This has to be worked out 
for each case separately and here we only give a couple of examples for illustration. 

Firstly and most qualitatively, assume that we have an SU(2)xU(l) theory without 
fermions. Using tree-level relations between physics and couplings and integrating out 
A Q , gives a simple analytic relation (eqs. (2.8-10) of 0]). According to that relation, 
x = 0.06444, 0.62402 correspond to m* H = 60, 180 GeV and z = 0, 0.3 correspond to 
m z = m w ,m e z P . 

Secondly, take the minimal standard model and the 1-loop equations relating x, y, z 
to physical parameters as given in Then x = 0.06444 corresponds to m# = 51.2 
GeV and x = 0.62402 to m H = 174 GeV. The value of z = z c at T c for given m#, m top 
is given in Fig. £| Since z c depends only logarithmically on T, we have estimated T c 
from y = (the true y c is very close to zero, see below). 

As a third example, take the minimal supersymmetric standard model. Fig. |^ shows 
different sets of values of the lightest Higgs mass uih, the right-handed stop mass m^ R 
and the CP-odd Higgs mass which all lead to the same value x = 0.06444. 

4 The lattice-continuum relations 

The lattice action corresponding to the continuum theory ( j2.1| ) is 

X i<j 
X i<j 
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Figure 2: The value of z = g'£/gi = z c at T = T c as a function of the physical 
Higgs mass m# and the top quark mass m top for the minimal standard model. The 
dependence shown arises mainly from the dependence of mw on m#, m top (see JlJ] for 
the formulas used), and gives an estimate of the uncertainty in the physical value of 
z c . The value of T c is computed from m\{g1) = 0; the dependence of z c on T c is only 
logarithmic. 



(4.1) 



Here the SU(2) and U(l) plaquettes are 



P v (x) 

Pij(x) 



Ui (x) Uj (x + i)Uj(x + j) U] (x) , 

exp{i[ctj(x) + oij{x + i) — ati(x + j) — aij(x)]}. 



(4.2) 
(4.3) 



Any positive number 7 gives the same naive continuum limit, but if exp(— 20^(2;) 03), 
plj" 1 are to be representations of the U(l) group and a G (0,27r), then one has to 
choose 7 = 1, 1/2, 1/3, . . .. The topological effects associated with 7 7^ 1 vanish in the 
continuum limit. 

A crucially important relation is the one expressing the four dimensionless lattice 
parameters /3a, l3' g , 0h, Pr of eq. ( |4.1[ ) in terms of the four dimensionless continuum 



parameters g%a, x, y, z. This relation, which can be derived using the techniques in |T6 
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x=0.06444 




180 200 220 240 

m~ 

t R 

Figure 3: Examples of parameter values corresponding to x — 0.06444 in the MSSM. 
Here mi R is the right-handed stop mass, inn is the lightest CP-even Higgs mass and 
rriA is the CP-odd Higgs mass. The squark mixing parameters have been put to zero. 
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1 - 2(3 R - 3P H ) + ^(3 + 12x + z 



32tt 



16tt 2 



51 



— - - ^z 2 + 9x - 12x 2 + 3xz^j (in ^(3 G + 0.09 



+5.0 - 0.9z + (0.01 + — )z 2 + 5.2x + 1.7x2 



(4.4) 
(4.5) 
(4.6) 



(4.7) 



where S = 3.17591. This relation is exact in the limit a — > 0. In [I7| it was further 
proposed (for 7 = 1) that the 1-loop 0(et)-corrections could be removed by modify- 
ing (|4.4|) - (|4~6D with terms of relative magnitude 0(1/ (3g)- 

In usual discretisations of the U(l) gauge theory, 7 = 1. In the physical situation the 
parameter z (eq.(|2.6|)) is about 0.3 and then the lattice coupling j3' G = j3 G /z (eq.( |4.6| )) 
would be rather large when approaching the continuum limit (3q — > 00. Using a 
value 7 < 1 permits one to have a smaller value of j3' G which makes finite size effects 
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smaller; the price one pays is an increasing correction term in eq. (|4.7|) and larger 0(a)- 
corrections related to g' 3 2 . 

The historical motivation for the form (|4.1|) of the lattice action was to permit one 
to go to the fixed length limit of the Higgs field by taking /3r — > oo. This limit is 
irrelevant now and it is actually possible very simply to combine (|4.1|) and (4~4])— ( |4.7|) 
by rescaling the Higgs field so that the coefficient of the quartic term is +1. The lattice 
action for fixed a and fixed continuum variables g 3 , x, y, z then is 



x i<j Z 



x i<j 



+ ^EE[i-i(ft/ 7 +4 /7 )] 



+ 



x 

Pg 

X 




8 S(3 + 12x + ^) 



1 



E[^Tr$t(a:)$(x)] 2 , 



E^Tr$ f (x)<l>(a;) 



(4- 



where the square bracket is the same as in ( |4.7| ). 

Finally, in the continuum limit, the relation of the gauge- invariant lattice observable 
(^Tr in ( J4.1| ) to the renormalized gauge-invariant continuum quantity {(jy(p(jj)) in 
the MS scheme is 



9l 



8tt 



167T 



^3 + z) (log ^ + 0.67 



2/i 



(4.9) 



The 1-loop 0(a) corrections to the discontinuity of (0V(aO) were also computed in |L7 



5 Perturbative results 



For understanding the magnitude of the expected effect of the U(l) subgroup, it is 
useful to compute it in perturbation theory. We have done this for the quantities listed 
in Section 2: the critical curve y = y c (x,z), the jump A£ 3 of the order parameter 
£ 3 = (<jy ~4>(gz)) I ' g\ between the broken (b) and symmetric (s) phases at y c and the 
interface tension a 3 . 

The coarsest approximation, which nevertheless gives the general pattern and can 
be given in analytic form, is obtained by taking the 1-loop potential and including only 
the vector loops. The result then is 
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Figure 4: The curves show the values of y c (x, z) computed from the 3d 2-loop effective 
potential of the SU(2) xU(l)+Higgs model. Lattice Monte Carlo results are also shown. 
The perturbative curves continue to large values of x, but on the lattice the first order 
transition ends at x ~ 1/8. Beyond that in the cross-over region, the values obtained for 
z = with the max(xi?2)-method (see Sec. 6.1) are y c (0.274, 0) = -0.066, y c (0.624,0) = 
-0.16. 



A4 



0"3 



128tt 2 x 2 
6(16vr) 3 UJ 



2\ 5 /2r2 1 



3 + 3 (1 + 2)3/2 



l 3 



(5.1) 



Here one sees concretely how increasing z increases the quantities studied. 

A more accurate result is obtained with 2-loop optimised 3d perturbation theory. 
The effective potential for the 3d SU(2)+Higgs theory is given and its optimisation is 
discussed in JTBJ. Since the effects of are very small, it is sufficient for the present 
purpose to differentiate between m 2 ^ and m 2 z only at 1-loop level (the complete 2-loop 
potential can be inferred from the 4d results in |19|, |2"0fl). A numerical computation then 
leads to the results shown in Figs. [|, |5|. The figures also give the results from lattice 
Monte Carlo simulations H || and make it clear how the perturbative discussion is a 
good qualitative quide, but fails in a crucial property of the transition: its termination. 
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Figure 5: The same as Fig. [|but for the jump of the order parameter and the interface 
tension. At x ~ 0.096, the interface tension was estimated in M to be 03 ~ 1 x 10~ 3 . 
For x > 1/8, there is no transition and A£ 3 , <t 3 vanish. 



6 Simulations 

The simulations with the SU(2) xU(l)+Higgs action ( |4.1|) are rather similar to those 
for the SU(2)+Higgs theory, described in detail in ||. We shall, therefore, mainly 
just present the results and compare them. We have performed simulations with x = 
0.06444 and x = 0.62402. For the SU(2)+Higgs case the former value gives a clear first 
order transition, whereas for the latter value the transition becomes a regular cross-over 
H £|. Since the most precise results for the SU(2)+Higgs theory were for x = 0.06444 
we shall concentrate mostly on this x-value. For z we always take z = 0.3. 

The simulations were mainly performed on a Cray C-90 supercomputer at the Center 
for Scientific Computing, Finland. Some of the smaller volume simulations were done 
on workstation clusters. The total amount of computing power used for the z = 0.3 
case was about 4.5 Cray cpu-months, or 2.3 x 10 15 floating point operations, to be 
compared with 5 x 10 15 flop for z — 0. 



6.1 The transition point y = y c (x = 0.06444, z = 0.3) 

Let us first discuss the x = 0.06444 case, where the transition is relatively strongly 
first order (m* H = 60 GeV in the notation used in ||). The SU(2)xU(l)+Higgs model 
does not have a local gauge invariant order parameter. We identify the transition by 
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studying order parameter like quantities, R 2 = J2x R 2 (x)/Vo\ and 

L ^ Wl ? ^ Tr ^WCaW* + Oe"***^ , (6-1) 

where is the SU(2)-direction of the Higgs field: $ = RV. These operators develop a 
discontinuity at the transition point, as shown by the probability distributions p{R 2 ) 
in Fig. |6| for (3 a = 5 and 8. The development of a strong two-peak structure with an 
increasing volume is unambiguous. 

In order to obtain the physical value of y c (x, z) we have to treat the finite volume 
effects and finite lattice spacing effects systematically. Our procedure here consists of 
the following stages: 

The pseudocritical coupling. For fixed x and (3q (lattice spacing) and each 
lattice size we determine the pseudotransition coupling (3 Hfi by several methods: 

(1) maximum location of the insusceptibility Xb? — ((R 2 ~ (R 2 )) 2 ), 

(2) minimum location of the Binder cumulant Bi = 1 — (L 4 )/(3(L 2 ) 2 ), 

(3) "equal weight" value for the distribution p(R 2 ), 

(4) "equal height" (3jj- value for p(L). 

In practice, we perform a series of simulations with different values of flu for fixed /3a 
and x (implying not fixed (3r but given by eq. (|4.5|)) until we have a good coverage 
of the transition region. The (3h- values corresponding to the above criteria are then 
found with the Ferrenberg-Swendsen multihistogram reweighting |HJ, and the error 



analysis is performed with the jackknife method. For all but the smallest volumes it is 
necessary to use multicanonical simulations For technical details we refer to [|3|, 
where a similar analysis was carried out for the SU(2)+Higgs model. 



X 


Pa 


7 


volumes 


0.06444 


5 


1 


8 2 x 32 m 10 2 x 40 m 
12 2 x 42 m 14 2 x 56 m 


8 


1 


24^ 32^ 16 2 x 112 m 20 2 x 120 m 
24 2 x 120 m 28 2 x 112 m 


16 


1/2 


12 3 16 a 24 a 32 a 

24 2 x 72 m 32 2 x 96 m 40 2 x 120 m 


0.62402 


8 


1 


12 3 16 3 24 3 32 3 



Table 1: The lattice sizes used for the simulations at the transition temperature for each 
(x, /3cf)-pair. In most of the cases, several /3^-values were used around the transition 
point. Multicanonical simulations are marked with the subscript ( m ). 

The lattice volumes used in the simulations are listed in Table |I|. For x = 0.06444 
we used /3q values 5, 8 and 16; in the last case we used 7 = 1/2 in order to avoid too 
large finite size effects. 
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SU(2)xU(l) x = 0.06444 z = 0.3 p G = 5 
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V y28 2 xll2 
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<7?> 

Figure 6: The "equal weight" distributions of the average Higgs length squared R 2 = 
Yj X R 2 {x)/V at the transition point for (5q = 5 (upper) and 8 (lower) for x = 0.06444. 

The infinite volume limit. For any given lattice, the criteria (l)-(4) above yield 
different values for the pseudocritical coupling (3h, c - However, in the thermodynamic 
limit V — > oo, all the methods extrapolate very accurately to the same value, as shown 
in Fig. [7[ It should be noted that the different methods do not yield statistically 
independent results, and combining the results together is not justified. In the extra- 
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SU(2)xU(l) x = 0.06444 z = 0.3 |3 G = 16 y= 1/2 

0.34072 ■ i ■ , — 

• max %(R ) 

■ equal weight R t 
o min B L A 
0.34071 □ equal height L I:: 



0.34070 



0.34069 O 



0.34068 



0.34067 1 1 1 1 1 1 1 1 1 

0.00000 0.00002 0.00004 0.00006 0.00008 

(volume/a ) 1 

Figure 7: The V = oo limit of the pseudocritical couplings (3h, c f° r x — 0.06444 
{m* H = 60 GeV) and 0a = 16, calculated with four different methods. The figures for 
Pg = 5, 8 are similar. 



Pa 


7 




Vc 


5 


1 


0.359863(3) 


0.00199(7) 


8 


1 


0.3489001(11) 


0.00255(7) 


16 


1/2 


0.3406899(8) 


0.00489(20) 


oo 




1/3 


0.0072(5) 



Table 2: The infinite volume critical couplings (3h, c an d the associated values of y c at 
x = 0.06444, z = 0.3 (from eq. (|4.7|) ). The values of @h, c are calculated from the "equal 
weight of p{R 2 ) n data. 



polations we used only volumes which are large enough in order to be compatible 
with the linear behaviour in 1/V. In Table [| we give the results using the "equal 
weight" values of the p{R 2 ) distributions. Also shown are the corresponding values of 
y c , obtained through eq. ( f4.7|) . Note that the value 7 = 1/2 for (3g = 16 enters only 
through a 2-loop contribution in eq. ( [4.7]) ; nevertheless, its inclusion is numerically 
crucial. 

The continuum limit. Lastly, the physical value of y c is obtained by extrapolating 
the y c (x,z;(3a) values to zero lattice spacing (p G = 4/(g%a) — > 00). This is shown in 
Fig. p. The leading behaviour is linear in a; however, it is clearly impossible to fit a 
straight line through all /?G-values. In principle, the finite a curve is different for each 
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Figure 8: The continuum limit (/3g — » cxo) of the critical temperature for x = 0.06444 
(m* H = 60 GeV), 2 = 0.3. 



value of 7, and fitting an extrapolation through (3g = 5 and 8 (7 = 1) and (3q = 16 
(7 = 1/2) datapoints is not rigorously justified. Nevertheless, when we compare this 
curve with the SU(2)+Higgs one (Fig. 11 in M), the qualitative similarity is obvious. 
However, the curvature of the quadratic fit is here clearly larger, and at (3q = 5 
the quadratic fit already starts to turn upwards. We assume that this behaviour is 
not physical and that the quadratic extrapolation overestimates the continuum value 
of y c . Therefore, in Table § we quote the value extrapolated linearly from (3q = 8 
and (3g = 16 -points; however, the difference between the linearly and quadratically 
extrapolated values is small and does not affect our conclusions. 

In Table |3] we compare the Monte Carlo results with the perturbative values. It 
is seen that for A£ 3 and cr 3 (whose determination is discussed below) the effect of z 
is just what is expected from the dominant 1-loop perturbative effect of g 3 2 . For y c 
the discrepancy is slightly larger; however, here the lattice determination of y c is not 
that easy (as explained in the previous paragraph). In any case the difference is small, 
and in terms of the critical temperature T c the discrepancy is below 1% since y c is 
composed of m 2 H and g 2 T 2 as a difference of two large terms. Thus for any physical 
conclusions, we can say that the last block in Table 10 of 0, which was based on a 
purely perturbative estimate of the effects of the U(l) group, can be considered reliable. 
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measured 


z = 


z = 0.3 


expected if U(l) perturbative 


tJc 


-0.00142(36) 


0.00724(45) 


0.0060 


. 2-loop 
i>C 


0.01141 


0.01882 






0.471(8) 


0.569(17) 


0.55 


A£ 2-loop 


0.493 


0.575 




_latt 
°3 


0.0116(28) 


0.0165(30) 


0.014 


2-loop 
^3 


0.0401 


0.0487 





Table 3: Comparison of measured values and expected values at x = 0.06444, z = 0.3 
if the effect of U(l) is perturbative. The 'expected' values are obtained by computing 
the difference (y c ) or ratio (A£ 3 , cr 3 ) of the lattice and perturbative values at z = 0, 
and by using these numbers in modifying the perturbative values for z = 0.3. Note 
that the z = values are slightly different from those in Table 10 of ||, since there the 
extrapolation to the continuum limit was made using the 4d observables T*, L/T* 4 , 
<j/T* 3 whereas here we use directly the 3d observables. 



6.2 The interface tension and A((ft(f)) 



We measure the interface tension with the histogram method ||23|| . At the transition 
point the probability of the two bulk phases is equal, but a system of a finite volume 
can also exist in a mixed state consisting of domains of the two bulk phases. The 
probability of the mixed state is suppressed by the free energy of the interface between 
the bulk phases. As a result, an order parameter distribution develops a characteristic 
2-peak structure (Fig. §), and the interface tension 03 can be extracted from the limit 



1 p max 

cr 3 = lim — ^ max (6.2) 



where P max an d -P m in are the values of the distribution at the peak and at the plateau 
between the peaks, and A is the smallest cross-sectional area on a three-dimensional 
periodic volume. In practice a careful finite size analysis is required; for technical 
details we again refer to the discussion in ||. 

The extrapolation of 03 to the infinite volume for each lattice spacing is shown in 
the upper part of Fig. |9| and to the continuum limit a = in the bottom part. The 
final result is shown in Table ^| together with the perturbative estimates. 

Let us next consider the measurements of the discontinuity of the condensate (<p'(/)). 
The lattice quantities R 2 = |Tr$^$ are first converted to the values of the continuum 
condensate at scale g\ by using eq. ( [4.9|) : £3 = ((p^ <p(g 2 )) jg\. The A{R 2 ) discontinuity 



is extracted from the "equal weight" distributions p{R 2 ) by measuring the positions of 
the peaks. The peak positions are determined by fitting parabolas to the distributions 
near the peaks. For each value of (3g the measurements are extrapolated to infinite 
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Figure 9: The interface tension cr 3 extrapolated to infinite volume (top) and then to 
the continuum limit (bottom) for x = 0.06444, z = 0.3. 



volume with respect to 1/A (Fig. |T0[) , and these values are in turn extrapolated to the 
continuum limit. The final results are shown in Table H. 
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Figure 10: The infinite volume limit of the discontinuity A(0^0)/g| for x = 0.06444, 
z = 0.3. The extrapolation is done linearly in 1/A, where A is the smallest cross- 
sectional area of the system. 



An alternative method for calculating A (</>V) is to use eq. ( p.l|) and measure 

A(0t0) -1 rfAlogZ -1 dAP 

9l Vg! dy Vgl dy ' ^ 

where AP is the difference of the probabilities of the two phases and the derivatives are 
evaluated at the transition point (where P symm . = -Pbrokcn = 1/2)- The result obtained 
with this method is quite compatible with the above one. 



6.3 Higgs and W correlators 

We measure the masses of H and W 3 with the lattice versions of the operators (jy(j) 
and Tr ^iD^a^, respectively. The correlation functions are measured in the direction 
of the 23 -axis, and in order to enhance the projection to the ground states, we use 
blocking in the (x\, X2)-plane. The fields are recursively mapped from blocking level 
(k) —►(& + !.) ($ (fc) (a;) — > $ (fc+1) (y)), so that the fields on the (k + l)-level lattice are 
in effect defined only on the even points of the (fc)-level lattice on the (x\, X2)-plane, 
doubling the lattice spacing along the (xi, x 2 )-plane: (ax, a 2 , a 3 )^ +1 ^ = (2a 1; 2a 2 , a 3 )( fc ). 
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Figure 11: m# and mw 3 measured from a. /3 G = 8, volume = 32 2 x 64 -lattice. Filled 
symbols correspond to broken phase, open symbols to symmetric phase masses. 



The blocking is performed with the transformations 

& k+1 \y) = I$( fc )( x ) + I ]T U$ k \x)$W(x + i)exp[-iai k \x)(T 3 \ (6.4) 

° ° t=±l,2 

and (i = 1, 2, j 7^ £) 

Ut l \y) = ut k \x)U« h \z + i), (6.5) 

^(x) = \ui%) + \ E ^(^^(x + jJ^Cx + i), (6.6) 
6 6 i=±i,a 

where (xi, X2, #3) = (2yi, 2y 2 , 2/3) and Z7_i(x) = U}{x — i). The U(l) gauge field S t - = 
exp[— zaf^03] is blocked with a transformation similar to ( |6.5| - |6.6| ). With the blocked 
fields we can define the operators 

h^(z) = E^Tr$«t( x )$W( x)? (6 .7) 

X 

wf k \z) = E^3,, Tr$Wt (^)^ (fc) ( a; ) $(fc) (^ + ^expHaf ) (a;) ( T3]^3- (6.8) 
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These operators are used to calculate the zero transverse momentum correlation func- 
tions 



= ^Y J h {k \z)h {k \z + l), (6.9) 
W^\l) = i-E w? k \z)w? k \z + l), (6.10) 

AV *;i=l,2 

and the masses are measured from the exponential fall-off in these functions. The 
blocking method used here resembles the one used by Philipsen et al. ||; however, 
we do not perform the diagonalization of the full cross- correlation matrix formed of 
operators of different blocking levels. 

For the Higgs channel the blocking makes little difference, and for the final results we 
use non-blocked (k = 0) operators. For W 3 , on the other hand, the blocking improves 
the signal substantially, and we use 3 times blocked operators W 3 ^ to measure the 
W 3 mass. At x = 0.06444 we measure the masses using /3q = 8, 32 2 x 64 lattices 
and performing simulations on both sides of the transition. The metastability of the 
transition is large enough to permit the mass measurements also in the metastable 
branches, as shown in Fig. 0. This figure can be compared with the SU(2)+Higgs 
masses in Fig. 16 of H; it should be noted that there the y range —0.066 < y < 0.042 



around the critical coupling y c is much wider than in Fig. 11 



6.4 The photon screening mass 

The qualitative difference between the SU(2)xU(l) theory considered in this paper 
and the SU(2)+Higgs model studied earlier is the presence of an extra gauge boson, 
associated with the U(l) group. This particle is massless in all orders of perturbation 
theory. 

There are several questions which can be addressed here. The first one is associated 
with the very existence of the massless vector excitation. Indeed, there are known 
examples when the photon gets a mass due to some non-perturbative effects. For 
example, in compact (discretized) QED in three dimensions the photon in fact does 
not exist and is replaced by a massive pseudo-scalar excitation EM. If some non- 



perturbative photon mass is generated, then, as has been argued in ||25|| , the problem 
of the primordial magnetic monopoles would be solved. 

The second question is associated with Z— 7 mixing. At sufficiently large and positive 
y (high temperatures) one would say that the massless excitation corresponds mainly 
to the hypercharge field B, while at large negative values of y (small temperatures) it 
is a mixture of the hypercharge field and the third component of the SU(2) gauge field. 
Since the SU(2)xU(l) theory does not contain a local gauge-invariant order parameter 
distinguishing the "broken" and "symmetric" phases, at sufficiently large values of the 
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parameter x (Higgs mass), the "content" of the hypercharge field in the massless state 
should interpolate smoothly between the two limiting cases. 

Both questions can be studied non-perturbatively on the lattice. In order to have a 
sensible plane-plane correlation function, we define the operator 



Op(z) = [Imui 2 (:Ei,22,z)]exp(zp ■ x), (6.11/ 



X 1 ,X 2 



where the sum is taken over the plane (xi, X2), Uij is the plaquette corresponding to the 
abelian U(l)-field, and p is a transverse momentum vector restricted to plane (1,2): 
ipiiP2iPz) — 2ir/N(ni, n,2, 0) with integer n« (there is quantization of the momentum 
because of the periodic boundary conditions). In the continuum limit the operator O 
is just 

O p (z) = / dxid,X2Bi 2 (xi, X2, z) exp(ip ■ x.) (6-12) 

^a(3' G J 

Note that if p = 0, eq. ( |6.12j) becomes identically zero in the absence of the winding 
modes. 

The correlator of two plane operators, 

G(z) = ^(YO p (t)0;(z + t)) (6.13) 



has the long-distance behaviour 



G{z) = A. QP exp(-^ v ^X), (6.14) 
z Pg Jv + mi 



where m 7 is the photon mass and Ay is measuring the projection of the lowest mass 
state to the hypercharge U(l) field. In the tree approximation we have v4 7 = 1 in the 
"symmetric" phase and = cos 2 ^ in the "broken" phase. For simplicity, eq. (|6.14|) 
is written in terms of the continuum dispersion relation instead of the lattice one; 
due to the small values of p used here the effect the effect of the different dispersion 
relations is well below the statistical errors in this analysis. As one can see, the use 
of a momentum dependent operator in Q6.ll ) is essential for a non-zero result for the 
correlation function. 

In the simulations we used the two smallest possible transverse momenta, p = 
2%/N ei, with i — 1 or 2. In Fig. [1^ we show the photon correlation function, measured 
from a x = 0.62402, (3g = 8 -system; all other parameter values yield very similar corre- 
lation functions. The results are perfectly consistent with the existence of the massless 
vector state. Applying eq. fl6.14p , the individual fits yield randomly positive or negative 
values for m 2 (compatible with zero within the statistical errors); as a final result we 
give an upper bound for the photon mass over the whole parameter range: 

m 7 < 0.034 (6-15) 
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Figure 12: The photon correlation function eq. ( 6.14Q for x = 0.62402 at the cross-over 
value for y, for a lattice of size 32 3 . The continuous line is a 2-parameter fit to the 
distance range 4-16. 



As to the mixing parameter Ay, in the vicinity of the phase transition with x = 
0.06444 we observe a sharp discontinuity in accordance with expectations (see Fig. |3~3j) . 
Note that here the discontinuity in is considerably smaller than given by the tree- 
level formula. This can to some extent be understood perturbatively; indeed, it is 
relatively easy to calculate at 1-loop order the residy of the 1/k 2 pole in the U(l) 
correlator (BiBj) (this determines Ay through (|6.12|) - (|6.14 )). In the symmetric phase 
(rajj > 0), one gets 

a 12 z 

A symm = 1 ^_ = 1 , g lg N 

7 487rm 3 48tt^/' V ; 

whereas in the broken phase the result is 

^broken = L + li^L\ ( 6 .1 7 ) 

' V 12nm w J 

where e\ = g 2 sin 2 ^^ = g\g'i / {g\ + g'i) and mw is the perturbative W mass. It can be 
verified that these expressions are gauge- independent. Both corrections have the right 
sign and even roughly the right magnitude. In the symmetric phase, the expansion 
breaks down close to the transition point as y ~ 0. 

As one approaches the endpoint of the line of the first order transitions, the disconti- 
nuity of A 1 gets smaller (together with my/). When x = 0.62402 (which is safely above 
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Figure 13: The pre-exponential factor A 7 for the x = 0.06444 case. The dashed 
horizontal lines are the tree-level values for in the symmetric and broken phases. 



the endpoint) we observe a continuous change in the mixing parameter Ay ~ 1 across 
the cross-over region. Thus, our non-perturbative analysis confirms within statistical 
errors the existence of the massless state which follows from perturbation theory, and 
demonstrates the absence of jumps of the mixing parameter in the region of large Higgs 
masses. 



6.5 Absence of spontaneous parity breaking 



A priori, it is not excluded that the Euclidean 3d non-abelian gauge-Higgs system 
exhibit spontaneous parity breaking in the strongly coupled phase. Arguments in favour 
of this possibility, based on the computation in 
topological mass term, were given in (28[ 



of the renormalization of the 3d 
Lattice simulations in the symmetric phase 
of the SU(2)+Higgs system were carried out in [29], where no signal for spontaneous 
parity breaking was found. Here we extend the analysis of |29| to the more realistic 
SU(2)xU(l) case. 

Consider some local gauge invariant operator O which is odd under the parity trans- 
formation. If there is spontaneous parity breaking in the system, the probability dis- 
tribution of the quantity 

^-fd 3 xO(x), (6.18) 

V •J 
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where the integral is taken over a finite volume V, must have a double-peak structure, 
which becomes more pronounced when the volume of the system increases. Simulta- 
neously, the susceptibility defined as 

X (V,T) = Jd 3 x(O(x)O(0)), (6-19) 

must behave as x(V, T c ) oc V 1 when V — > oo. If spontaneous parity breaking is absent, 
the probability distribution of O looks like a single peak with center at zero, and 
X (V,T c )<xV°. 

Our choice of parity breaking operators includes the the lattice versions of the oper- 
ators ( [2.15| - f2~T6| ) and the operator 



W P (x) = e ijk W t 3 (x)W+(x)W l :(x) 1 (6.20) 



where W ± = W 1 ± iW 2 . The study of the non-local operators (|2.17|j correspond- 



ing to the Chern-Simons numbers, was not attempted because it would be very time 
consuming. 

We found that the expectation values of all of these operators remain consistent 
with zero in both symmetric and broken phases at x — 0.06444 near the phase tran- 
sition. More precisely, the distribution of the operators is very accurately Gaussian 
centered around zero; moreover, the autocorrelations of the successive measurements 
are negligible even in the transition region, where all of the other measurements show 
significant autocorrelations. This implies that the parity operators decouple from the 
long-distance dynamics and are sensitive only to the local ultraviolet noise of the sys- 
tem. The susceptibility \ was found to be practically volume independent for all of 
the volumes, again signifying the non-critical nature of the parity operators. The sit- 
uation is similar at x — 0.62402 in the cross-over region. Thus, the possibility of 
non-perturbative parity breaking is ruled out (at least within the statistical accuracy 
and for the volumes used in our simulations) for the SU(2)xU(l) theory as well. 

6.6 The transition in the x = 0.62402 case 

When x = 0.62402 (in this value of x corresponds to m* H = 180 GeV) and z = 0.3 
the transition has turned into a smooth and regular cross-over, completely analogously 
to the SU(2)+Higgs z = case. This can be observed, for example, by monitoring 
the behaviour of the order parameter susceptibility and the correlation lengths. For 
example, in Fig. [T4I we show the behaviour of the ftcj) susceptibility \<t> as a function 



of y for several volumes; in the case of a true transition the susceptibility diverges at 
the transition point (oc V for a first order transition, cx V 1 for a second order one). In 
this case no divergence is seen, only a volume independent smooth peak corresponding 
to the cross-over behaviour. For this x we have data only for (3q = 8, and using the 
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Figure 14: The </>'(/) susceptibility x<t> measured from (3q = 8, x = 0.62402 lattices 
as a function of y. No divergence develops when the volume increases (top). The 
maximum values of the susceptibility function of the volume for SU(2)+Higgs 

and SU(2)xU(l)+Higgs systems (bottom). 



maximum location of x<f> we extrapolate the infinite volume cross-over coupling to be 
y c (x = 0.62402, z = 0.3, (3 G = 8) = -0.152(7). 

In Fig. [14] we compare the maximum values of x<i> f° r SU(2)+Higgs and SU(2) xll(l) + 
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Figure 15: The scalar correlation length from a (3 a = 8, 
lattice as a function of y. The dotted line shows the location of the maximum of x<j>- 



Higgs theories at x — 0.62402; the SU(2)+Higgs data is from J4|. The quantitative 
behaviour of the systems is quite similar, and there is no first order transition. In 



Fig. |T| we show the (lightest) scalar mass across the transition region; it remains non- 
zero as in the SU(2)+Higgs, signalling the absence of a second order transition. The 
photon mass vanishes everywhere in both phases (see Sec. fT3), and hence it should 
not affect the order of the transition. 



7 Conclusions 

We have in this paper studied numerically the 3d SU(2) xU(l) -(-fundamental Higgs 
theory and compared it with the 3d SU(2) -(-fundamental Higgs theory at the same 
values of the coupling x = A 3 /g|, but for ^ 0. The comparison has been made 
both in a region with a clear first-order transition (x = 0.06444; m* H = 60 GeV) and 
in the cross-over region (x = 0.62402; m* H = 180 GeV). The main conclusion is that 
the inclusion of U(l) does not bring in significant new non-perturbative effects: the 
changes in the non-perturbatively computed quantities due to increasing g'^ are as large 
as those given by perturbation theory. At the same time, the values of the observables 
themselves, especially 03, differ from those in perturbation theory. In particular, the 
first order transition ends at some x c as in 3d SU(2)+Higgs and for x > x c there is 
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no first or second order transition. A still higher order transition cannot be strictly 
excluded by the analysis presented in this paper. 

We have also confirmed non-perturbatively the existence of a (nearly) massless pho- 
ton in the 3d SU(2) xU(l)+Higgs model. We have studied the mixing of this state to the 
Lagrangian U(l)-field, and demonstrated that the jump in the mixing parameter van- 
ishes together with the first order transition as one goes to large enough Higgs masses. 
Physically these results on the photon imply that magnetic fields are not screened in 
hot electroweak plasma even taking into account all nonperturbative effects. 
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